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Методом теории случайных процессов установлены вероятности пребывания системы в различных 
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цессе эксплуатации, что, в свою очередь, обеспечит надежность эксплуатации технологического обо- 
рудования в различных отраслях промышленности. 
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Постановка проблемы 


Уровень развития современной науки и техники существенно зависит от внедрения 
новых материалов. Поэтому особенно актуальны исследования свойств эпоксидных 
композитных материалов (ЭКМ), которые являются лучшими заменителями тради- 
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ционных материалов во многих отраслях промышленности [1], [2]. Широкое исполь- 
зование ЭКМ стало причиной интенсивного исследования их физико-механических 
свойств и разработки методов их прогнозирования. Важнейшим фактором повышения 
эффективности производства в любой отрасли является улучшение управления про- 
цессами сшивки материалов и их функционирования под воздействием внешних ста- 
тических и динамических напряжений. 

Функционирование широкого класса систем можно представить как процесс пе- 
рехода из одного состояния в другое под воздействием каких-либо причин. Применим 
один из эффективных математических методов теории случайных процессов, в част- 
ности, один из ее разделов — теория марковских процессов [3], [4] для решения задачи 
прогнозирования свойств полимерных композитов. 


Анализ последних исследований и публикаций 


При исследовании непрерывных и дискретных случайных цепей обычно поль- 
зуются графическим представлением функционирования системы. Граф состояний 
системы представляет собой совокупность вершин, изображающих возможные состоя- 
ния системы 5», и совокупность ветвей, изображающих возможные переходы системы 
из одного состояния в другое. Марковский случайный дискретный процесс, протека- 
ющий в системе 5, характеризуется не только возможными состояниями, в которых 
система может пребывать случайным образом, но и теми моментами времени, в ко- 
торые могут происходить ее переходы из состояния в состояние. Эти моменты вре- 
мени могут быть заранее известны или случайны [5]. 

Случайный процесс, протекающий в системе, называется процессом с дискретным 
временем, если переходы системы из одного состояния в другое могут осуществляться 
только в заранее определенные моменты времени 11, №,......, называемые шагами 
этого процесса. В промежутках между соседними шагами система сохраняет свои 
состояния. Не исключается возможность, что на некоторых шагах система не изменит 
своего состояния. 

Случайная последовательность называется марковской цепью, если для каждого 
шага вероятность перехода из любого состояния 5, в любое состояние 5; не зависит 
от того, когда и как система 5 оказалась в состоянии 5. 

Так как система 5 в любой момент / может пребывать только в одном из состоя- 


ний 5) и, ТО при каждом К = 1.2,... события К (^).. р КУ (К) несовместны и обра- 


зуют полную группу. 
Цель исследования — определить вероятности состояний системы (эпоксидного 
композитного материала) при воздействии статической нагрузки. 


Постановка задачи 


На предварительном этапе в результате эксперимента получена зависимость абсо- 
лютной деформации от продолжительности воздействия статической нагрузки эпоксид- 
ного композита, содержащего 50 масс.ч. наполнителя (электрокорунд) на 100 масс-ч. 
эпоксидной смолы. 

При исследовании ползучести ЭКМ использовали стандартную методику испыта- 
ний на изгиб по ГОСТ 4648-71 при экспериментально установленной статической 
нагрузке Г = 50 Н. Схема устройства для исследования ползучести материалов при- 


356 «Искусственный интеллект» 2013 №3 


Определение вероятностей напряженного состояния... 4к_ 


ведена на рис. 1. Отметим, что исследовали материалы в течение времени т = 72 ч, 
при этом определяли абсолютное значение прогиба образца под нагрузкой в конкре- 
тно определенные моменты времени. Параметры образцов: длина — [= 120 мм, ши- 
рина - р = 15 мм, высота — й = 10 мм. 


Рисунок 1 — Схема устройства для исследования ползучести материалов: 
1 — несущая нагрузка; 2 — индикатор; 3 — идентор; 4 — образец; 5 — опора 


Поставим задачу: найти вероятности состояний системы. Для решения этой за- 
дачи построим математическую модель исследуемого процесса. 

Изменение абсолютной деформация (ЛГ) под воздействием статической нагруз- 
ки ЭКМ во времени (1) позволяет установить следующую серию состояний: Зо— исходное 
состояние; 5; — работоспособное состояние; 5 — упругая деформация; 5з;— пластичная 
деформация; 5 — разрушение (табл. 1). Марковский случайный дискретный процесс, 
протекающий в системе 5, характеризуется не только возможными состояниями, в ко- 
торых система может пребывать случайным образом, но и теми моментами времени, 
в которые могут происходить ее переходы из состояния в состояние. Эти моменты вре- 
мени могут быть заранее известны или случайны. 


Таблица | — Зависимость абсолютной деформации от времени воздействия 
нагрузки ЭКМ 


Состояния системы 
50 $1 52 53 54 


Ч 0,0112 4 6 12 | 18 | 24 | 30 | 36 42 48 | 54 | 60 | 66 | 72 
4[х10°, м| 68 | 71 | 74 | 76 | 79 | 80 | 81 | 83 | 84 84,6 | 85 86 87 |88 |90 


Параметр 


Очевидно, поставленная задача относится именно к марковской цепи, поэтому 
граф состояний может быть представлен в виде (рис. 2). 


“— 5—5 ——5. 
Юл 


Рисунок 2 — Граф состояний системы 
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Основными характеристиками марковских цепей являются вероятности р, (Ё)= 
= ($. (®))@=1,.. пу =1,2,...) событий 5, (К). Вероятности р, (К) (1=1...., т; К =1.2....) 
являются вероятностями состояний. Таким образом, вероятность 71 состояния на К шаге 
Р.(Ё) является вероятностью того, что система 5 от А до (&+/) шага будет пребывать 


в состоянии 5;. Сумма вероятностей этих событий для каждого К = 1,2,... равна 1: 


Ур... 
1=1 


Пусть система 5 в момент / находится в состоянии 5,. Рассмотрим элементарный 
промежуток времени 1, примыкающий к моменту 1. Назовем плотностью вероятности 
перехода 1, из состояния 1 в состояние / (или инфинитезимальными коэффициентами) 
предел отношения вероятности перехода системы за время 1 из состояния 5; в со- 
стояние 5; к длине промежутка /#: 


(+ -р.( ‚(А 
хора +А0- РО РУ 
1—0 МЕ 4—0 УАУА 


ый 


где Р‚(Д - вероятность того, что система, находившаяся в момент / в состоянии 5, 
за время /1 перейдет из него в состояние 5, (плотность вероятностей перехода опре- 
деляется только для 7). Отсюда следует, что при малом 4 вероятность перехода (с точ- 
ностью до бесконечно малых высших порядков) равна: Р‚(И) =, 1. 

Если все плотности вероятностей перехода 1, не зависят от { (от того, в какой 
момент начинается элементарный участок 1), марковский процесс называется одно- 
родным, а если эти плотности зависят от времени, то он является неоднородным. 

Запишем вероятности состояний в соответствии с графом состояний (рис. 2) в 
виде квадратной матрицы и порядка, сумма элементов каждой строки равна 1 (рис. 3). 

Дискретный случайный процесс с дискретным временем, протекающий в системе, 
характеризуется тем, что система может переходить из одного состояния в другое толь- 
ко в заранее определенные моменты времени, называемые шагами. 

Известно, что множество состояний системы марковской цепи определенным обра- 
зом классифицируются с учетом дальнейшего поведения системы следующим обра- 
зом [6]. 


Рисунок 3 — Теоретическая матрица состояний 


1 Невозвратное множество. 

В случае невозвратного множества возможны любые переходы внутри этого 
множества. Система может покинуть это множество, но не может вернуться в него. 

2 Возвратное множество. 
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В этом случае также возможны любые переходы внутри множества. Система 
может войти в это множество, но не может покинуть его. 

3 Эргодическое множество. 

В случае эргодического множества возможны любые переходы внутри мно- 
жества, но исключены переходы из множества и в него. 

4 Поглощающее множество. 

При попадании системы в это множество процесс заканчивается. 

Кроме описанной выше классификации множеств различают состояния системы: 

а) существенное состояние: возможны переходы из 5; в 5; и обратно; 

6) несущественное состояние: возможен переход из 5, в 5, но невозможен обрат- 
НЫЙ. 

В некоторых случаях, несмотря на случайность процесса, имеется возможность 
до определенной степени управлять законами распределения или параметрами пере- 
ходных вероятностей. Такие марковские цепи называются управляемыми [6]. 

Кроме того, с учетом наличия или отсутствия тех или иных упомянутых выше 
множеств состояний марковские цепи могут быть поглощающими, если имеется хо- 
тя бы одно поглощающее состояние, или эргодическими, если переходные вероятности 
образуют эргодическое множество. 

Практически важным является вопрос о том, сколько шагов сможет пройти система 
до остановки процесса, то есть поглощения в том или ином состоянии. 

Очевидно, согласно построенному графу состояний системы (1), состояния 50, 
51, 52 являются существенными состояниями, так как возможны переходы из 5, в 5, 
и обратно (1.1 = 0,1,2); состояние 53 является невозвратным множеством, так как систе- 
ма может покинуть это множество, но не может в него вернуться. Состояние 54 является 
поглощающим множеством. 

Вычислим вероятности состояний Р,. Поскольку в данном случае воспользоваться 
классической формулой вероятности Р = и/и невозможно, определим вероятности 
состояний по частоте наблюдения того или иного значения наблюдаемого признака. 

Определение вероятности по частоте наиболее широко применяется в науке 
и технике, когда не выполняются условия для непосредственного подсчета по форму- 
ле Р = иии. Частоту определяют по статистическим данным или по результатам спе- 
циально проведенных опытов. В нашем случае находим число т (время наблюдения 
конкретного состояния), и общее число и (общее время наблюдения за процессом в 
тех же условиях); тогда вероятность состояния Р* = пии. 

По табл. 1 вычислим время 1; пребывания системы в каждом конкретном состоянии 
абсолютной деформации ЛГ, (табл. 2). 


Таблица 2 — Время пребывания системы в различных состояниях при воздействии 
статической нагрузки ЭКМ 


Состояния системы 
п 
мет р |5 52 53 54 
-3 
ый 68 [71| 74 | 76 | 79 | 80 | 81 | 83 | 84| 846| 85 | 86 | 87 | 88 | 90 
ат м И И" 


Если рассматривать весь процесс в целом, общее время наблюдения 
(суммарная частота) составляет (= 72 ч. 
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Обозначим (1; время пребывания системы в состоянии 1, соответственно абсолют- 
ная деформация состояния Л/;. 
Тогда вероятности состояний Р; абсолютной деформации АГ, за весь период наблю- 


м 
хе 


Р езультаты вычисления вероятностей состояний приведены в табл. 3. 


дения равны р, = 


Таблица 3 — Вероятности пребывания системы в состояниях 50... 54 


Состояния Параметры 
системы АЁ;х1 0 м Ч Частота, А Ёч Вероятность 
50 68 0,1 0,1 0,0014 
$1 71 2 1,9 0,02639 
74 4 и. 0,0278 
76 6 р. 0,0278 
79 12 6 0,0833 
г 80 18 6 0,0833 
81 24 6 0,0833 
83 30 6 0,0833 
84 36 6 0,0833 
84,6 42 6 0,0833 
53 85 48 6 0,0833 
86 54 6 0,0833 
87 60 6 0,0833 
88 66 6 0,0833 
54 90 12 6 0,0833 


Поскольку процесс марковский, будем рассматривать отдельные этапы перехода 
процесса из одного состояния в другое, т.е. каждый этап будет представлять смежные 
состояния: 50, 51, 52 (первый этап, восстанавливаемая деформация); 5/, 52, 53 (второй 
этап, переход в пластичную деформацию); 53, 54 (третий этап, невосстанавливаемая 
деформация). Три смежных состояния (50, 51, 52) наблюдаем в течение {= 24 ч; деформа- 
цию А, наблюдаем с вероятностью 0,04167; деформацию Д/> - с вероятностью 0,07917 
и т.д. (табл. 4). Так как эти события независимы (в каждый момент времени система 
находится только в одном состоянии), вероятность каждого из состояний вычисляем 
как сумму вероятностей составляющих этого состояния (табл. 4). 


Таблица 4 — Вероятности пребывания системы в состояниях 50... 52 


Параметры 
Состояния _3 Сумма вероят- 
АЁ,х10”, м 
системы ы Ч Частота, ч Вероятность | ностей в 
состояниях 
50 68 0,1 0,1 0.04167 0,04167 
71 2 1,9 0,07917 
ы! 74 4 2 0,0833 и 
76 6 2 0,0833 
79 12 6 0,25 
2. 80 18 6 0,25 м 
81 24 6 0,25 
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График зависимости абсолютной деформации от продолжительности воздействия 
статической нагрузки представлен на рис. 4 


-3 Состояния $0, $1, $2 
АИ `,м 


[=] 
ю 


80 
78 
76 


74 


абсолютная деформация 


72 
70 


58 
56 т 


время наблюдения 


Рисунок 4 — График зависимости абсолютной деформации от продолжительности 
воздействия статической нагрузки 
Аналогично вычислим вероятности состояний третьей и четвертой строк матрицы 
состояний (табл. 5, табл. 6) и построим графики зависимости деформации (рис. 5, рис. 6). 


Таблица 5 — Вероятности пребывания системы в состояниях 5... 53 


Параметры 
а Сумма вероят- 
системы А И х] 0" М Е ч Частота, р Версия ее 
СОСТОЯНИЯХ 
71 2 1,9 0,02883 
— 74 4 2 0.030349 0,05918 
76 6 2 0,03035 
79 12 6 0,091047 
-- 80 18 6 0.09105 0,30349 
81 24 6 0,09105 
83 30 6 0,0909 
84 36 6 0,0909 
84.6 42 6 0.0909 
53 85 48 6 0,0909 0,6364 
86 54 6 0.0909 
87 60 6 0,0909 
88 66 6 0,0909 
АТО м Состояния $1 ‚ $2. $3 
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Рисунок 5 — График зависимости абсолютной деформации от продолжительности 
воздействия статической нагрузки 
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Таблица 6 — Вероятности пребывания системы в состояниях 53... 54 


Параметры 
Состояния Сумма вероят- 
системы АГ, х 103, м Ч Частота, ч Вероятность ностей в 
состояниях 

83 30 6 0,125 

84 36 6 0,125 

84,6 42 6 0,125 

53 85 48 6 0,125 0,875 

86 54 6 0,125 

87 60 6 0,125 

88 66 6 0,125 

54 90 72 6 0,125 0,125 
3 С остояния $3.54 
АТА ^,м 
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Рисунок 6 — График зависимости абсолютной деформации 
от продолжительности воздействия статической нагрузки 


Таким образом, запишем матрицу состояний (табл. 7). 


Таблица 7 — Матрица состояний системы при воздействии 
статической нагрузки ЭКМ 


50 51 52 53 54 Сумма 
50 0 1 0 0 0 1 
51 0,04167 0,1625 0,8333333 0 0 1 
52 0 0,05918 0,30349 0,63733 0 1 
53 0 0 0 0,875 0,125 й, 
54 0 0 0 0 1 1 


Для получения более полной характеристики состояний системы вычислим ма- 
тематические ожидания деформации под воздействием статической нагрузки ЭКМ, а 
также дисперсию, т.е. разброс относительно наиболее ожидаемых значений Деформа- 
ции. Как известно [6], математическим ожиданием (средним значением) случайной 
величины Х, заданной на дискретном вероятностном пространстве, называется число 
т = М/Х] = Ухр,. Результаты вычисления математического ожидания состояний каждого 
этапа марковской цепи приведены в табл. 8. 

Таблица 8 — Математическое ожидание состояний марковской цепи 


Состояния системы Математическое ожидание, М(Х) Дисперсия, О(Х) 
50... 54 83.43472х10 м 12,40907х103 м 

50... 52 78.40417х10 м 9.29915х103 м 

53 $2.86039х103 м 15,9587х10 м 

54 $5,95000х10° м 4,6175х10° м 
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Определение вероятностей напряженного состояния... 4к 


Полученные значения имеют следующий смысл. При воздействии статической 
нагрузки К = 50 Н в композитах, содержащих 50 масс.ч. наполнителя электрокорунда 
на 100 масс.ч. эпоксидной смолы на первом этапе 50... .52 в течение / = 24 ч наиболее 
вероятна абсолютная деформация М(Х) = 78, 410 мс дисперсией Д(Х) = 9,3х107 г 
при воздействии этой же статической нагрузки в течение 1=60ч наиболее вероятна 
абсолютная деформация М(Х) = 82, 9х103 мс дисперсией П(Х) = 15, 96х10 ; при не 
ствии от 30 до 72 часов наиболее вероятна абсолютная деформация м = 86х103 
с дисперсией Д(Х) = 4,6х 10° 
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Деаегтттаноп о{РгоБаб Иез 5тез5 5кие оГЕроху Сотро5йе 


Маетав ипаег “айс Гоаа 

ш изше фе ШФеогу оЁ МаКоу ргосеззез а шефо4 г рге@сйпе сВапоез ш фе 
ргорегие$ оЁ ро[утег сотрозЦез \Пеп ехрозе4 1ю $айс 1оа4. 

АЕ Фе ргейпитагу %азе оЁ Фе ехрегитепе ргодисе а рйом шЮюптайоп абойе 1е 
свапое оГ шаена| ргорегиез аз а РапсНоп оР абзоцие Аеогтайоп оЁ Фе диганоп оЁ 
ехрозиге фо Зайс |оа4 оГероху сотрозпе. Апа[у$1$ оЁ Фе гези5 аПоууе4 11$ 10 еза $6 
Фе зресблит оЁ за ез оРШе зу%ет апа сопзвтасЕ а огарь оРа @1зсгее МагКоу гапдот ргосезз. 

Озше Фе Шеогу оЁ зюсвазИс ргосеззез ргофабИиу оЁ Ше зубет Беше шзаПе ш 
уапопз за ез, сопзгасе тафлх оЁ йе зузет зез ипдег фе шИчепсе оЁ Фе $айс |оа4 оп 
ерок$Котро7йпу таепа! НПе ул рагяез Разед. Тпе 4еуе]оре4 тоде| аПоууз ю ргедс 
Фе ргорегаез ап Бепау1ог оЁ сотрозИе таепа[5 ип4ег 5ез$ Чиние Фе ореганоп, УЛисв, шт 
бит, епзиге тепаЫе орегайоп оЁ ргосезз едиртеп ш уапои$ шаизлез. 
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